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We discuss applications of the quantile concept of trajecto- 
ries and velocities to the propagation of electromagnetic sig- 
nals in wave guides of varying cross section. Quantile motion 
is a general description of the transport properties of measur- 
able conserved quantities in quantum mechanics as well as in 
classical electrodynamics. In the latter case we consider the 
quantile motion of electromagnetic energy as the direct result 
of a physical measurement. In this sense the quantile veloc- 
ity corresponds to the electromagnetic signal velocity also in 
the presence of barriers and inhomogeneities in the medium 
of propagation. We show that this signal velocity is always 
smaller than the speed of light in vacuum. Using numerical 
examples we demonstrate how typical wave phenomena can 
be described in terms of the quantile motion. 



I. INTRODUCTION 

In recent years the tunneling of an evanescent electro- 
magnetic pulse in a wave guide with a cross section re- 
duced over part of the longitudinal extension of the wave 
guide has been studied. A series of microwave experi- 
ments claims to have found superluminal tunneling 
velocities. The velocity of the signal has been identified 
with the velocity of the maximum of the pulse. 

A physical detector always needs a certain amount of 
energy to change its state, therefore we use the quantile 
velocity of a density obtained by way of normalizing the 
electromagnetic energy density of the pulse to unity || |t| . 

A numerical simulation of the microwave experiments 
with wave packets in wave guides with varying cross sec- 
tion is carried through. The quantile velocity of wave 
packets is studied numerically. It is shown rigorously 
that the quantile velocity is not superluminal. A short 
discussion of the quantile velocity in dispersive and ab- 
sorptive media is presented. For a critical discussion of 
superluminal velocities of electromagnetic signals see also 
P. Thoma at al. [§, G. Raithel j|, G. Diener |Tf| and 
the contribution of H. Goenner in this workshop 



II. STATIONARY WAVES IN WAVE GUIDES OF 
VARYING CROSS SECTIONS. UNITARITY 
RELATIONS 

We shall consider a wave guide of rectangular cross 
section and perfectly conducting walls, extending in the 
longitudinal z direction from — oo to +oo. The simplest 
construction allowing an interpretation in terms of tun- 
neling for an electromagnetic signal propagating in the 
wave guide is obtained by assuming a narrowing of the 
cross section, e.g., in the transverse x direction extending 
to the left of the origin for z > as shown in Fig. VI. 
An electromagnetic wave moving from the right towards 
the narrowing of the cross section will be affected in a 
similar way as the quantum-mechanical wave function of 
a particle incident onto a potential step at z = 0. 

To further simplify the problem we consider only the 
propagation of TE- waves of the type H„o for n = 1, 2, . . .. 
Such a situation can indeed be realized if we choose a con- 



stant height of the wave guide (b in Fig. VI) small enough 
for all waves of the type H nm with m > 1 to be evanes- 
cent Thus, the stationary electromagnetic field in 
the wave guide reduces to three field components, E y , 
H x , and H z , each depending on x and z. One of the field 
components can be considered as independent, e.g., E y , 
whereas the other two follow from the Maxwell equations, 
e.g., iw/i/ioH = V x E, where ui is the frequency of the 
wave and /i = 1 is the magnetic permeability of vacuum. 
It is worth noting that according to a well-known result 
p3| , p^| the magnetic field components are singular at the 
sharp edge of the narrowing. This implies an increase of 
the magnetic energy at the cost of the electric energy in 
the wave guide [|l2|,[l4| ■ 

The solution for the stationary electric field component 
E y fulfilling the Helmholtz wave equation with the appro- 
priate boundary conditions is a superposition of modes 
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where we denote k„ = y^(uj/c) 2 — (nir/a) 2 and n' n — 
yj{u)/c) 2 — (nir/a') 2 the discrete wave numbers associ- 
ated with the two regions of the wave guide. Considering 
a given mode H m o a s the incoming wave incident from 
the right (z < 0) we are led to the ansatz for the solution 
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for z < and z > 0, respectively. Here, .A^ and A^ are 
the (complex) amplitudes of the reflected and transmit- 
ted modes in complete analogy to the one-dimensional 
tunneling in quantum mechanics. However, an infinite 
number of different modes is needed in order to fulfill the 
boundary conditions at z = 0, so that the type of the 
outgoing wave does not necessarily coincide with that of 
the incoming wave. 

Requiring continuity for the electric field component 
E y at the wave-guide aperture at z = 0, we obtain for 
the amplitudes A^ and A^ 
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where n = 1,2,... and e(x) = E y (x, z = 0) is the value 
of the electric field strength at the aperture. Denoting 
S mn +A n - = A£> and A„ = Aff we can equivalently write 
equations ((H) as 
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where 
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The continuity condition for the transverse magnetic 
field component at z — yields 
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Introducing equations (||) into (|5|) we obtain for the elec- 
tric field e(x) at the aperture the integral equation 
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for < x < a. The kernel IC(x, x') is given by 

= E»=i [*? sin(^) sin(^) 
+ 4sin(2^)sin(^) . 
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The reflection and transmission amplitudes .A^ and 
follow then from the equations (||) and the solution of the 
integral equation ([]]). 



In general the integral equation (g) is not solvable in 
closed form and the field e(x) has to be determined nu- 
merically. We can, however, try to find approximate ex- 
pressions for A„ and A„ by using the alternative repre- 
sentation of the integral equation (Q) as a system of linear 
equations for the coefficients A^ , 
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with k = 1, 2, . . . and matrix elements T n fc given by 

oo 
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Under certain conditions discussed below we can assume 
the matrix T to be nearly diagonal and put its elements 
in the form 



Tnk = «4( y ) 2 (#nfc + m n fe) , 



(10) 



where m is a complex matrix with elements fulfilling 
\mnk\ <C 1 for n, fc = 1, 2, . . .. The solution of m) follows 
then from the expansion of the corresponding inverse ma- 
trix into a fast converging Neumann series 
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Neglecting all the terms in ( |11|) but the leading one we 
obtain the coefficients 
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for n = 1, 2, . . ., implying at z 
field component H x given by 



a transverse magnetic 
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Equation (|lj) predicts the increase in the magnitude of 
the transverse magnetic field component H x at z = 0, 
whereas the shape of the field there obviously coincides 
with that of the incoming mode H.,„o . Thus, the approx- 
imation leading to (12) seems appropriate in the limiting 
case of geometric optics, i.e., if the wave length of the 
incoming wave is short compared to the cross sections a 
and a' of the wave guide. In the case of an evanescent 
wave for z > 0, however, more terms in the expansion 
( |TT| ) are needed for the boundary conditions at z = to 
be fulfilled. Therefore, a numerical solution of equation 
(^) is used in what follows. 

We now turn to a more complicated geometry for the 
wave guide involving two boundary conditions at z — 
and z — L, where L is the length of a symmetrically 
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placed barrier (reduction of the cross section as shown in 
Fig. VI). Numerical simulations for the stationary field 
components in the wave guide are shown in Fig. VI and 
Fig. |]. The continuity of the transverse field components 
at the wave-guide aperture as well as the magnetic-field 
singularities at the narrowing edges are shown in Fig. [vj|. 
Figure |^ is a spectral diagram with reflection and trans- 
mission coefficients (in the regions z < and z > L, 
respectively) of the incoming as well as of higher modes. 

Because of the magnetic-field singularities in a real 
experiment the maxima of the magnetic field strength 
would always be located in the vicinity of the wave-guide 
edges. Therefore, they give no information on the veloc- 
ity of a tunneling electromagnetic signal. If, instead, we 
refer to the electromagnetic energy density in the wave 
guide we can avoid dealing with singularities by recalling 
that the singularities of the energy density are integrable. 
Poynting's theorem in our case reads 
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where u> om is the two-dimensional energy density and s x 
and s z are the components of the Poynting vector. In- 
tegrating dl^) over the transverse direction x the second 
term in (|14|) vanishes and we obtain 
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where W(z,t) and S z {z,t) are now the one-dimensional 
energy and current densities in the longitudinal direc- 
tion z. These quantities are free of singularities and 
fulfill the one-dimensional continuity equation (^|) in 
complete analogy to the probability and current density 
in the one-dimensional tunneling in quantum mechanics 
[pTsfj . Therefore, it is appropriate to consider the energy 
density W(z,t) instead of the electric and magnetic field 
strengths in order to investigate the tunneling properties 
of electromagnetic signals in wave guides. 

Returning to the stationary fields, we may use jl5| ) 
to obtain unitarity relations between the reflection and 
transmission coefficients of the modes A„ and n = 
1,2,.... We find the time- averaged longitudinal current 
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to be a constant along the wave guide. Comparing the 
expressions for S z {z) in the reflection and transmission 
region we obtain 
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where H m o is the incoming mode, H„ c o and H„/o are the 
cutoff modes in the reflection and transmission region, 
respectively. The appearance in ( p"7j ) of the upper limits 



n c and n' c in the summation index n is a consequence of 
time averaging the current density S z {z). Thus, equation 
( |i~7| ) does not imply that only propagating modes n < n c 
(n < n' c ) are responsible for the energy transport through 
the barrier. In the microwave experiments |l|-^| the cross 
sections a and a' of the reflection and transmission region 
were the same. With this assumption equation (J17|) reads 
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Further simplifying the situation by allowing the ground 
mode Hio as the only propagating mode in the wave 
guide yields the unitarity relation of the form 
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which clearly coincides with the corresponding equation 
in the one-dimensional quantum-mechanical tunneling. 



III. TUNNELING OF WAVE PACKETS. 
COMPARISON WITH EXPERIMENTS 

In the following we construct tunneling wave packets 
in the wave guide as a superposition of the stationary 
solutions found above. For such a superposition to be a 
propagating, incoming wave in the region to the right of 
the barrier (z < 0) it must not contain wave components 
below the cutoff frequency lu c (evanescent components) 
in this region. Therefore we use, e.g., for the electric field 
component E y the expression 
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where f(uj) is the spectral distribution of frequencies nor- 
malized in the interval (u c , oo) and E ys is the stationary 
electric field component. 

In Fig. ||-Fig. |?j numerical simulations with Gaussian- 
like wave packets in the wave guide are shown. The spec- 
tral distribution f{uo) in this case is given by 
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where N is a normalization factor, loq is the mean fre- 
quency and (Ju] is the spectral width of the wave packet. 

Figure ^ demonstrates the tunneling process through 
a symmetrically placed barrier in the wave guide for an 
incoming mode of the type Hio- Most of the frequencies 
in the spectrum ( ^l[ ) are taken above the cutoff frequency 
uj' c of the barrier region L > z > 0. Correspondingly, the 
transmission rate is high. We observe that the shape of 
the reflected and transmitted wave packets is substan- 
tially deformed. Among other reasons this is due to the 
energy transfer to higher modes other than the incoming 
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mode Hio- In terms of the evolution of the wave-packet 
maxima, we start with a single maximum in the incoming 
wave and end up with many transmitted maxima propa- 
gating with different velocities in the wave guide. 

In Fig. ^| the resonant tunneling of a Gaussian wave 
packet is shown. The frequency spectrum contains many 
transmission resonances (see also Fig. ^) leading to mul- 
tiple reflections of the wave packet at the barrier walls. 

In Fig. ^ we compute the longitudinal energy density 
W(z, t) in the time domain at a fixed position z behind 
the barrier and for different barrier lengths L. The wave 
packet was chosen such that tunneling takes place mainly 
in the evanescent regime. Under this condition we ob- 
serve that the maximum of the wave packet can appear 
behind the barrier earlier than when moving in the free 
space with the vacuum speed of light c. This behavior 
becomes more obvious for large L as the transmission 
rate decreases significantly. 

As another example we consider the tunneling of 
Kaiser-Bessel wave packets [|l6) with a limited and dis- 
crete frequency spectrum in a given interval [cj_,c<j + ]. 
The spectral distribution is given by 
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where a is the parameter characterizing the width of the 
distribution and N is the number of the stationary-wave 
components with frequencies uj n = lu- + (n + -y)(o;+ — 
lu-)/N for n = —N/2, . . . , —I, 0,1..., N/2. The distri- 
bution (^2|) ensures optimality of the wave-packet local- 
ization in the time domain and was also used in the 
microwave experiments [Q. In Fig. ^ we compute the 
tunneling time of the maximum of a Kaiser-Bessel wave 
packet between the beginning and end of the barrier as a 
function of the barrier length L. We consider wave pack- 
ets tunneling in the propagating and evanescent regime 
and compare in the respective case the tunneling ve- 
locity of the maximum with the group velocity in the 
wave guide and the velocity of light in vacuum c. In 
the evanescent case we observe the tunneling velocity of 
the maximum of the wave packet for long barriers to be 
independent of the barrier length L. This behavior cor- 
responds to the Hartman effect which is well-known from 
the quantum-mechanical tunneling |l7],[ll|. Thus, evanes- 
cent tunneling of the maximum becomes highly superlu- 
minal. In the example of Fig. [s] we obtain a tunneling 
time tt ~ 127.78 ps at a barrier length L = 100 mm 
for the maximum of the wave packet. The correspond- 
ing tunneling velocity is «t = L/tt « 2.6c in very good 
agreement with the experimental result given in Q| . 



IV. QUANTILE MOTION OF 
ELECTROMAGNETIC SIGNALS. CAUSALITY 
OF SIGNAL PROPAGATION 

We measure the arrival time of a signal with a detec- 
tor placed at the fixed position z P . We assume that the 
region in which the energy of an electromagnetic pulse is 
essentially different from zero is initially far away from 
the position z P of the detector. The detection of the elec- 
tromagnetic signal requires the deposition of a certain 
amount W of energy in the detector to cause a change 
of its state indicating the arrival of the signal. This is 
equivalent to the condition 



dzW{z,t) = W 
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on the time t of arrival of the signal. Repeated mea- 
surements at different positions Zpi, z P 2, • ■ ■ yield arrival 
times t\,t2,.. - corresponding to these positions. They 
are discrete points on the trajectory z P = z P (t) defined 
by requiring the condition 
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dzW(z,t) = W 
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to hold at all times t. If we call Wq the total energy 
contained in the pulse then P = W/Wo is the fraction of 
energy needed for detection and g(z,t) — W(z 7 t)/Wo is 
the normalized energy density. Equation ( 24 ) can be put 
into the form 



zp(t) 



dzg(z,t)=P, 0<P<1 
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which is the same as Equation (11) of |J. Therefore, 
z P = z P (t) is the quantile trajectory of the electromag- 
netic signal. As to be expected it depends on the fraction 
P, and thus on the sensitivity W — PWo of the detector. 
The signal velocity is then given by 



>(*) = 



dz P (t) 
dt 
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Examples of quantile trajectories for tunneling Gaus- 
sian signals in wave guides of varying cross section are 
given m Fig. [v| and Fig. [vj. Figure ^ shows that the 
presence of a barrier in the wave guide may only lead 
to a slower signal propagation at any give time t and 
for every detector sensitivity P. Thus, no quantile signal 
velocity larger than the speed of light in vacuum c is pos- 
sible. Especially, in the evanescent tunneling regime (see 
Fig. |Vl|c) the tunneling velocity is much smaller than c, 
whereas, as physically expected, most of the trajectories 
turn back to the reflection region. 

The behavior of the quantile trajectories for different P 
values reflects several properties of the tunneling process. 
In Fig. |V4^ 3 the reflected and transmitted quantile tra- 
jectories split into trajectory bunches propagating with 
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different velocities in the wave guide. They correspond 
to the electromagnetic modes produced in the tunneling 
as described in the previous section. This coincidence 
between the behavior of quantile trajectories and typi- 
cal wave phenomena in tunneling can be observed also in 



resonant tunneling (see Fig. VI). 

The causality of quantile motion can explicitly be de- 
rived, for instance, in the case of tunneling TE-waves. 
Using the definition ( |24| ) and the continuity relation ( |lE| ) 
we derive the quantile velocity fl24|), c.f. Equation (13) 
in §, 



v P (t) 



dz q (t) S[z P (t),t] 



dt 
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Trajectories solving this equation have also been stud- 
ied by Holland |19| . The modulus of the velocity field 
characterizing the differential equation ( p7|) is 



\v(z,t)\ 



=1 Re[cE« 



Re[H2] 



< c. 
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Thus, v q (t) never exceeds the vacuum speed of light c, 
i.e., the signal propagation described by the quantile tra- 
jectory is causal. This result is a general property of the 
quantile motion and holds independently of the type of 
the tunneling wave. 



exchange is characterized by a loss of energy, since part of 
the mechanical energy of the oscillators is steadily trans- 
formed into heat. Thus, three different kinds of energy 
have to be considered, the electromagnetic energy of the 
pulse, the mechanical energy of the oscillating charges 
of the medium, and the thermal energy stored in the 
medium p7| ]. Each of them may be considered as the 
energy of an open subsystem of the closed system char- 
acterized by the total energy which is conserved. 

We apply the concept of quantile motion to the wave 
propagation in dispersive media. Sufficient condition for 
this is the existence of a measurable quantity with a pos- 
itive spatial density (see (^-0]). Since the total energy in 
the medium is conserved, its density fulfills a continuity 
equation similar to (|l4|). Thus, quantile trajectories for 
the total energy in a dispersive and absorptive medium 
may be easily defined in complete analogy to the above 
description for the (conserved) electromagnetic energy in 
wave guides. The application of quantile motion for the 
electromagnetic signal propagation in dispersive media 
has been carried through in |27j] . In all cases it has been 
shown that the quantile velocities for the signal remain 
below the vacuum speed of light c. An argument anal- 
ogous to the one leading to (^8|) shows rigorously that 
the quantile velocity in the oscillator model of dispersive 
and absorptive media is always smaller than the vacuum 
speed of light. This and related results will be presented 
in a forthcoming publication. 



V. NOTE ON THE QUANTILE MOTION IN 
DISPERSIVE AND ABSORPTIVE MEDIA 



It has been known for a long time that electromag- 
netic signal propagation in the spectral region of a disper- 
sive medium characterized by anomalous dispersion and 
strong absorption leads to superluminal phase and group 
velocity [p0|-p3|. Even though the shape of the propa- 
gating signal may be substantially deformed in compar- 
ison to the shape of the incoming wave, a propagation 
velocity has been considered which coincides with the 
velocity of one of the pulse maxima. Then the velocity 
obtained this way is again superluminal p3] | . This result 
was recently reconfirmed by the photon experiments at 
Berkeley |24|,^j| and by the experiments with photonic 
barriers in wave guides p6| . However, if one considers 
instead of the pulse maxima the energy transport in the 
medium no superluminal velocities occur. Diener fic| l 
suggests a procedure for separating propagating and non- 
propagating parts of energy of the electromagnetic pulse 
and obtains subluminal velocities for the energy trans- 
port in the medium. 

Because of the interaction of the electromagnetic field 
with the dispersive medium, energy is permanently trans- 
ferred between the propagating pulse and the oscillating 
charges in the medium. In the presence of dissipation this 



VI. CONCLUDING REMARKS 

The concept of quantile motion has been applied to 
the propagation of electromagnetic waves in wave guides 
of varying cross section and in dispersive and absorptive 
media. It has been shown that the signal velocity mea- 
sured with a detector of finite sensitivity never becomes 
superluminal. In the context of the propagation of elec- 
tromagnetic signals the quantile velocity is a generaliza- 
tion to detectors with finite sensitivity of Sommerfeld's 
concept M] of a front velocity describing the speed of a 
signal measured with a detector of infinite sensitivity. 
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FIG. 1. Geometry of a wave guide with a reduction of the 
cross section beginning at z = and extending to infinity. 
This configuration with an H„o wave incident from the right 
corresponds to a quantum-mechanical potential step. 




FIG. 2. Geometry of a wave guide with a symmetric reduc- 
tion of the cross section between z = and z = L correspond- 
ing to a barrier of finite length L. 




FIG. 3. Stationary fields in evanescent tunneling. The di- 
mensions of the symmetric wave guide and of the barrier are 
a = 1 cm, a = 0.4 cm, and L = 1 cm (see Fig. |vi| ). The incom- 
ing wave incident from the right is an Hio mode of frequency 
u — 21 GHz. (a) The electric field strength \E y ^~ is shown. 
This component is continuous everywhere and vanishes at the 
walls of the wave guide, (b) The transverse magnetic field 
strength l-H^] 2 is shown. Obviously, the magnetic field is sin- 
gular at the edges at z = and z = L of the wave guide. 
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FIG. 4. Spectral diagram with the transmission and reflection coefficients of the first first two modes Hio and H20, where 
H10 is the incident mode. The dimensions of the wave guide and of the symmetric barrier are a = 1 cm, a' — 0.5 cm, and 
L = 1cm (see Fig. Q). 




FIG. 5. Time development of a tunneling Gaussian wave packet. The electric field strength \E y \ 2 is shown. The dimensions 
of the wave guide are a = 1cm, a — 0.5 cm, and L = 5 cm. The incident wave packet is of the type Hio with frequencies 
centered at coo = 30 GHz and spectral width — 1.5 GHz. The production of higher modes in tunneling which propagate with 
different velocities in the wave guide, is observed. Times are given in picoseconds. 
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(a) t = (b) t = 11.67 (c) t = 23.33 




(d) t = 35.00 (e) t = 46.67 (f) t = 58.33 




FIG. 6. Time development of a tunneling Gaussian wave packet in resonant tunneling. The electric field strength \E y \ 2 is 
shown. The dimensions of the wave guide are the same as in Fig. The spectral function of the incoming wave packet centered 
at u>o = 18.97 GHz with a u = 0.015 GHz extends over more than one resonance (see Fig. so that multiple reflection inside 
the barrier is observed. Times are given in picoseconds. 




FIG. 7. Example of the superluminal tunneling of the maximum of a tunneling Gaussian wave packet. The longitudinal 
energy density W(z, t) is shown in the time domain at a fixed position behind the barrier Az — 20 cm away from the position 
of the center of the wave packet at t = 0. The dimensions of the symmetric wave guide are a = 1 cm and a = 0.5 cm for 
different barrier lengths L (see Fig. VI). The incident wave packet is of the type Hxo with frequencies centered at u>o = 15 GHz 
and spectral width = 0.6 GHz. Thus, evanescent tunneling occurs. The box, the triangle, and the circle correspond to the 
arrival times of the pulse maximum in vacuum, in a wave guide without barrier, and in the wave guide with a barrier of length 
L, respectively. If the barrier is long enough the tunneling velocity of the pulse maximum becomes superluminal (case (c) and 
(d)). 
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FIG. 8. Tunneling times of the pulse maximum and Hartman effect. The transmission time of the pulse maximum (solid 
line), the time intervals corresponding to the free motion of the pulse maximum in vacuum (dashed line) and in the wave 
guide without barrier (broken line) for a Kaiser-Bessel wave packet are plotted as functions of the barrier length L. In both 
cases the incoming wave packet is a superposition of N = 800 frequencies in the interval 51.52 GHz < uo < 57.81 GHz with the 
spectral distribution ( p2] ) for a = 1. (a) Non-evanescent tunneling. The cross sections of the wave guide are a = 22.86 mm and 
a' — 18 mm. The behavior of the tunneling times is "normal" . (b) Evanescent tunneling. Here, the same tunneling situation 
as in is computed, i.e., the two cross sections of the wave guide are a — 22.86 mm and a' = 15.8 mm. The transmission time 
r T = L/v m for long barriers remains nearly constant and, thus, implies superluminal tunneling for the pulse maximum. 
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FIG. 9. Quantile trajectories, (a) Non-evanescent tunnel- 
ing. Trajectories for 0.05 < P < 0.95 in steps of AP = 0.02 
are shown. The dimensions of the symmetric wave guide (see 
Fig. VI) are a = 1 cm, a' = 0.3 cm, and L = 5 cm. The incom- 



ing Gaussian wave packet is of the type Hio with u>o = 60 GHz 
and Ou) = 0.3 GHz. (b) Evanescent tunneling with short bar- 
rier. The cross sections a and a 1 are the same as in (a). 
The length of the barrier is L — 0.1cm. The incoming wave 
packet with luq = 30 GHz and = 0.3 GHz is evanescent 
in the region inside the barrier, (c) Evanescent tunneling 
with long barrier. The length of the barrier is chosen to 
L — 1 cm. The incoming wave packet has u>q = 22.5 GHz 
and (Ju, — 0.3 GHz. Because of the strong reflection, trajecto- 
ries for 10" 3 < P < 2.5 x 10~ 2 in steps of AP = 5 x 10~ 4 
are computed. Also in this case no superluminal tunneling is 
observed. 
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FIG. 10. Quantile trajectories in resonant tunneling. The 
dimensions of the symmetric wave guide are a = 1 cm, 
a' = 0.5 cm, and L = 100 cm. The incoming wave is a 
Gaussian wave packet of the type Hio with coo = 18.97 GHz 
and cju, — 0.015 GHz. (a) Time development of the longi- 
tudinal energy density W(z,t). (b) Quantile trajectories for 
0.2 < P < 0.625 in steps of AP = 0.025. The oscillating 
behavior of the quantile trajectories reflects the time develop- 
ment shown in (a). 
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